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Abstract The pressure on reducing environmental footprint is facilitating the 
emergence of energy supply chains that have biomass as main feedstock. For the 
development of efficient energy supply chains from biomass it is required to 
properly integrate the various elements that comprise such systems (e.g., biomass 
supply, pre-treatment facilities and technologies for biomass to energy and/or fuels 
conversion). Additionally, it is recognised that a concerted effort is required, 
embracing the different supply chain entities, in order to correctly estimate envi¬ 
ronmental burdens and to propose effective environmental strategies. This chapter 
proposes the use of a mathematical modelling approach as an analytical tool that 
can support decision-making towards accomplishing the design and planning of 
efficient multiple source—multiple product bio-energy supply chains. The math¬ 
ematical formulation of this problem becomes a multi-objective MILP (moMILP). 
Criteria selected for the objective function are the net present value (NPV) and the 
overall environmental impact, which is computed using the Impact 2002+ 
indicator. The main advantages of this approach are highlighted through a case 
study of a biomass-based supply chain that comprises different components 
geographically distributed over Spain. For comparison purposes, such a supply 
chain is contrasted to one embracing coal usage. 
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Notation 

Indices 

e Suppliers 

/, f Facility locations 

i Tasks 

j Equipment technology 

s Materials (states) 

t , t' Planning periods 

a Mid-point environmental impact categories 

g End-point environmental impact categories 

Sets 

A g Set of mid-point environmental interventions that are combined into 
end-point damage factors g 
E rm Set of suppliers e that provide raw materials 

£prod Set of suppliers e that provide production services 

E tr Set of suppliers e that provide transportation services 

E e Set of locations / where supplier e is placed 

FP Set of materials s that are final products 

Ij Set of tasks i that can be performed in technology j 

J e Technology j that is available at supplier e 
Jj Technology j that can be installed at location / 

J t Technologies that can perform task i 

Mkt Set of market locations 

NTr Set of non-distribution tasks 

RM Set of materials s that are raw materials 

Sup Set of supplier locations 

T s Set of tasks producing material s 

T s Set of tasks consuming material s 

Tr Set of distribution tasks 


Parameters 


^sft 

Dem sft 
distance^ 


FCFI 


j J 

L ft 


] jft 


Norm. E g 
Pricey 

Pric 4 


Maximum availability of raw material s in period t in location / 
Demand of product s at market / in period t 
Distance from location / to location / 

Fixed cost per unit of technology j capacity at location / in period t 
Investment required to establish a processing facility in location / in 
period t 

Normalising factor of damage category g 
Price of product s at market / in period t 

Investment required per unit of technology j capacity increased at 
facility / in period t 
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rate 

Water, 

Water*/* 

Msij 

&sij 

h 

U 

Ow 


tiff't 


T 


utl 

ijfet 


r ut2 

T sfet 


tf]a 


Discount rate 
Moisture for material s 

Maximum moisture for task i performed in equipment j 
Mass fraction of task i for production of material s in equipment j 
Mass fraction of task i for consumption of material s in equipment j 
Minimum utilisation rate of technology j capacity that is allowed at 
location / 

g end-point damage characterisation factor for environmental inter¬ 
vention a 

Capacity utilisation rate of technology j by task i whose origin is 
location / and destination location f 

Unitary transportation costs from location/to location / during period t 
Unitary cost associated with task i performed in equipment j from 
location / and payable to external supplier e during period t 
Unitary cost associated with handling the inventory of material s in 
location / and payable to external supplier e during period t 
Unitary cost of raw material s offered by external supplier e in period t 
a environmental category impact CF for task i performed using 
technology j receiving materials from node/and delivering it at node/ 
a environmental category impact CF for the transportation of a mass 
unit of material over a length unit 


Binary Variables 

Vjft 1 if technology j is installed at location / in period t , 0 otherwise 


Continuous Variables 


DamC^ 

DamC^ c 

EPurch^ 

ESales ? 

FAsseh 

FCosh 

Fjft 

FE jft 

aft 

Impact^ 002 

Impact 2 ^ 

LHV« 

NPV 

Pijfft 


Normalised end-point damage g for location / in period t 
Normalised end-point damage g along the whole SC 

Economic value of purchases executed in period t to supplier e 
Economic value of sales executed in period t 
Investment on fixed assets in period t 
Fixed cost in period t 

Total capacity of technology j during period t at location / 
Capacity increment of technology j at location / during period t 
Mid-point a environmental impact associated to site / which rises 
from activities in period t 
Total environmental impact for site/ 

Total environmental impact for the whole SC 
Lower heating value for material s in task i 
Net present value 

Activity magnitude of task i in equipment j in period t whose origin 
is location / and destination location / 
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Pv si jff t Amount of material s for flexible task i in equipment j in period 

t whose origin is location / and destination location f 

Profit^ Profit achieved in period t 

Purch^ r Amount of money payable to supplier e in period t associated with 

production activities 

Purch^ Amount of money payable to supplier e in period t associated with 
consumption of raw materials 

Purch^ Amount of money payable to supplier e in period t associated with 

consumption of transport services 

Sales 5 ^/ Amount of product s sold from location / in market f in period t 

Ssft Amount of stock of material s at location / in period t 


Superscripts 

L Lower bound 
U Upper bound 


Acronyms 

B-NET 

Biomass utilisation networks 

BM 

Biomass 

CBA 

Cost benefit analysis 

CCS 

Carbon capture and storage 

CF 

Characterisation factors 

DFCF 

Discounted-free-cash-flow 

Eco-indicator 99 

Damage environmental metric 

EISA 

Energy Independence and Security Act 

FWR 

Forest wood residues 

GHG 

Greenhouse gas 

GrSCM 

Green supply chain management 

IMPACT 2002+ 

Mid-point and end-point (damage) environmental metric 

LCA 

Life cycle assessment 

LCIA 

Life cycle impact assessment 

LCI 

Life cycle inventory 

STN 

State task network 

moMILP 

Multi-objective mixed integer linear programming 

NPV 

Net present Value 

SCM 

Supply chain management 

STN 

State task network 
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1 Introduction 


Global trends are promoting the utilisation of renewable sources of energy. 
For instance, the European Union has established a target of 20% share of 
renewable energy out of the total European energy consumption by 2020 [1]. 
Biomass energy generation is an option that is expected to play an important role 
in the renewable energy mix. As a matter of fact, the USA in its Energy 
Independence and Security Act (EISA) of 2007 states that advanced biofuels shall 
supply at least 21 billion gallons of US motor fuels by 2022. It is recognised that in 
order to achieve these targets efficient networks to supply sustainable amounts of 
biomass required, cost-effective technologies to convert biomass and improved 
distribution infrastructures to deliver the final product (i.e., energy or fuel) are to 
be developed [2]. Even more, the efficient integration of these three elements is 
equally relevant to reach the posed goals. In this chapter, we propose a supply 
chain modelling approach as a tool that can support decision-making towards 
accomplishing such integration. 

The concept of supply chain (SC) refers to the network of interdependent entities 
(i.e., processing sites, distributors, transporters, warehouses and raw material 
suppliers) which is the processing and distribution channels of a product from 
origin of its raw materials to final delivery to the customer. Then, supply chain 
management (SCM) can be defined as the management of material, information and 
financial flows through a SC that aims at producing and delivering goods or services 
to consumers [3]. Note that a SC is comprised by components that may be 
geographically distributed. One of the main objectives of SCM is to synchronise 
and coordinate the flows of materials that go through the different processes so that 
the final product is delivered in a most efficient manner. This is especially important 
for biomass to energy projects which are highly geographically dependent and 
whose profitability can be strongly influenced by the location of the different 
processes and biomass sources. Commonly, biomass production and transportation 
account for a significant part of the whole bioenergy supply chain cost [4]. 
Therefore, a tool capable of evaluating the possible trade-offs between the different 
feedstock sources, each one with specific properties (i.e., moisture and energy 
density) and the location of processing sites and consumption is a requisite to 
develop efficient bioenergy networks. 

The interest on SC biomass-based for different final purposes, ranging from 
energy or fuels production, has increased since the 90s, thus being a relatively new 
concept and also a consequence of the current energy paradigm. Typically, a 
biomass SC problem considers the possible use of multiple sources of biomass. It 
is also important to point out that biomass usually requires some pre-treatment 
before its use. This previous treatment is done aiming at obtaining a homogenised 
or densified biomass, whether in terms of mass or energy. In this sense, Hamelinck 
et al. [5] considers an international bioenergy SC taking into account the fact that 
biomass production and consumption do not need to be in the same region, i.e., 
that biomass residues and energy crops are in excess in some areas, but scarce in 
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others, and energy import regions (Western Europe) do not match with the higher 
biomass production areas (Scandinavia, Eastern Europe and Latin America). 
Biomass compacted into briquettes or pellets can be used to save transport use 
because of their smaller volume. Nevertheless, there exists a trade-off between the 
type of transport, the distance and the densification methods since the cost of 
biomass pre-treatments as a whole should not be higher than the compensated 
logistic costs. This is an important consideration to take into account when long 
distances should be covered. They demonstrate that an international bioenergy 
trade has real potential; however, governments must stimulate the biomass market 
in terms of prices, policies and social acceptance. Panichelli and Gnansounou [4] 
contemplate forest wood residues (FWR) from final cuttings to produce torrified 
wood that supplies a gasification unit in order to produce electricity. They are able 
to allocate biomass quantities between predefined combinations of candidate sites 
to find the best set of locations for the energy units. They state fix values for the 
torrefaction and gasification units, and take into account cost minimisation. A 
mixed integer linear program (MILP) that determines the optimal sizes and 
locations of biomass-based methanol plants (biofuel plants) is developed by Leduc 
et al. [6]. The objective function considered to be optimised is the operating costs 
plus the investment required to establish the biomass SC. The supply is given by 
poplar coppice, as energy crop, and the demand is based on gasoline-methanol car 
blend use. The possible consumer sites are the already existing gas stations in 
Austria, while methanol production is considered through the use of gasification 
plants. The evaluation of three scenarios is performed, based on different metha¬ 
nol-gasoline blends. By-production of heat is also considered as economic reve¬ 
nue, and C0 2 emissions are accounted finally, but not introduced in their model as 
an environmental objective to be accomplished. Later on, Ayoub et al. [7] present 
a methodology for designing and evaluating the biomass utilisation networks (so 
called B-NETs), which are process networks aiming at producing different bio¬ 
products, from one or more biomass resources. The idea of this methodology is to 
provide a framework to create the underlying superstructure that relates the bio¬ 
mass resources to their products via current and possible future available pro¬ 
cesses; which can be used to develop an optimisation model. Their methodology is 
applied at a local level and proposes better biomass uses by means of economic 
parameters (costs) and environmental impacts accounting. The environmental 
indicator that the authors use accounts for emissions to air, water pollutants and 
solid wastes. The work of Rentizelas et al. [8] emphasises the multi-biomass 
seasonal availability and combines this fact with the biomass storage problem. The 
stages considered before the conversion plant of the raw material include har¬ 
vesting and collection, in field handling and transport, storage, loading and 
unloading, transport and biomass pre-treatment. This last stage can be included in 
any of the abovementioned stages, and optimally could precede the transportation 
stage. Storage can be also located in the biomass origin in an intermediate step or 
at the power station site. The authors state that one of the main drawbacks of the 
use of biomass as a source is its relatively low density and heating value when 
compared, for instance, with other fossil sources. On a later work, Rentizelas et al. 
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[9] exemplify the fact that a biomass SC can account for multiple sources, as well 
as for multiple final products production, such as electricity, heat and cooling. 
They apply the methodology in a specific region of Greece. The results provide 
with optimal locations and investment details for potential investors. More 
recently, Van Dyken et al. [10] developed a linear optimisation model for planning 
the capacity expansions in energy systems where several alternative biomass and 
technologies are considered simultaneously. The main objective of this work is to 
present a generic model including different components such as sources, handling, 
processing, storage and final usage. Heating value, moisture content and bulk 
density are the key parameter changes that biomass undergo along the SC. The 
objectives to be optimised are the operating cost and emissions of the whole SC. 

Definitely, energy policies are driven by environmental considerations, more 
specifically by the pressure on reducing greenhouse gas (GHG) emissions. In that 
sense, biomass is an energy source that is expected to provide significant reductions 
of environmental impacts related to GHG emissions when compared to the classical 
fossil fuels technologies. Therefore, it is relevant the integration of environmental 
thinking into SCM in order to assess such expected reduced environmental impacts. 
The aforementioned integration may be achieved through the emerging concept 
regarded as ‘green supply chain management’ (GrSCM). This concept considers the 
environmental interventions associated with the raw materials sourcing and selec¬ 
tion, manufacturing process selection, delivery of final product to the consumers as 
well as end of life management of the product after its useful life [11]. Traditionally, 
the methodologies devised to assist SC operation and design have focused on finding 
a solution that maximises a given economic performance indicator while satisfying a 
set of operational constraints imposed by the manufacturing/processing technology 
and the topology of the network. In recent years, however, there has been a growing 
awareness of the importance of including environmental aspects as objectives and 
not constraints associated with the SC decision support [12, 13]. 

The environmental science and engineering community have developed several 
systematic methodologies for the detailed characterisation of the environmental 
impacts of chemicals, products and processes. All of these methodologies have 
embodied the concepts of life cycle, i.e., they are based on a life cycle assessment 
(LCA) which is described in a series of ISO documents [14]. The LCA framework 
includes the entire life cycle of the product, process or activity, encompassing 
extraction and processing of raw materials; manufacturing, transport and distri¬ 
bution; re-use, maintenance recycling and final disposal. Most importantly, it takes 
a holistic approach, bringing the environmental impacts into one consistent 
framework, wherever and whenever these impacts have occurred or will 
occur [15]. These methodologies are based on the incorporation of an optimisation 
step into the four classical phases that comprise an LCA study namely, goal 
definition, life cycle inventory—LCI, life cycle impact assessment—LCIA and 
interpretation (see Fig. 1). The idea of them is to determine process conditions or 
topology using a multi-criteria optimisation strategy in order to evaluate the trade¬ 
off between economic and environmental issues. 
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Fig. 1 Life cycle assessment 
steps [13] 
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As aforementioned, the concept of SC refers to the network of interdependent 
entities that constitute the processing and distribution channels of a product from 
the supply of its raw materials to its delivery to the final consumer. Because an 
LCA study ideally covers a cradle-to-grave approach, it can be clearly seen that 
LCA fits as a suitable tool for quantitatively assessing the environmental burdens 
associated with designing and operating a SC. 

This chapter describes an analytical approach for the design and planning of a 
multiple source—multiple product interregional bioenergy SC taking into con¬ 
sideration not only economic issues but also environmental impact. The model 
accounts for different biomass wastes sources as well. The approach applies mixed 
integer modelling techniques. The model is optimised so as to select the most 
appropriate pre-treatment technologies, the most appropriate feedstock supplier 
for each process and the most convenient production and distribution profiles, in 
order to supply electricity and hydrogen to the customers. The mathematical model 
encompasses direct emissions, raw materials production and transport distribution 
emissions. LCA concepts are embedded in the approach, and going further in order 
to attain a comprehensive LCA application, not merely an overall environmental 
impact indicator is calculated but also partial environmental impact categories are 
studied. Furthermore, the impact associated with every SC echelon is mapped 
aiming at discovering possible opportunities to focus resources for environmental 
impact reduction. 
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2 Biomass Supply and Operations Planning 

This chapter deals with the strategic-tactical problem associated to the optimal 
design and operation of a biomass SC network taking into account economical and 
environmental considerations. Generally speaking, the SC strategic level 
determines the network through which production and distribution serves the 
marketplace. The intent of the SC network design problem is typically to deter¬ 
mine the optimal sourcing, manufacturing and distribution network for the new 
and existing product lines of a company (e.g., expansion or contraction of the 
business, introduction of new products, new strategic suppliers). The most 
common approach is to formulate a large-scale mixed integer linear program 
(MILP) that captures the relevant fixed and variable operating costs for each 
facility and each major product [16]. 

The considered biomass SC network consists of a number of potential locations 
where either a processing site or distribution centre or both of them can be 
installed, and suppliers at fixed locations which have available biomass with 
different characteristics. In general, each final product (energy or H 2 ) can be 
produced at several plants located at different locations using the different biomass 
wastes. The characteristics of the biomass can be changed by using the 
pre-treatment units (e.g., drying or torrefaction) so that the treated biomass meets 
the characteristics required to be used in further steps. Even more, such pre¬ 
treatments increase the energy content and bulk density of the biomass. By doing 
so the mass is reduced, thus significant savings in transportation may be achieved. 
The production capacity of each processing site is modelled by relating the 
nominal production rate per activity to the availability of the equipment tech¬ 
nology at each plant. Distribution centres are described by upper and lower bounds 
on their material handling capacity and they can be supplied from more than one 
manufacturing plants. Given the way the problem is modelled, materials flow 
between any facilities may appear if selecting such flow allows improving the 
performance of the SC. A market demand may be served by more than one site. 

The mathematical model is an analytical tool intended to support managers on 
planning decisions such as: 

• The active SC nodes and links among them 

• The facilities capacity expansion in each time period 

• The product portfolio per plant, production amounts, utilisation level and 
transportation links to establish in the network alongside with material flows 

• The amount of final products (energy or H 2 ) to be sold 

• The environmental impact associated to each SC node or activity 

A general schematic of the biomass energy SC is shown in Fig. 2. One can 
notice that it is comprised by four blocks: (i) sourcing, (ii) pre-treatment, 
(iii) generation and (iv) distribution. The sourcing block consists in collecting the 
different biomass that may be available from different regions and suppliers. 
Each type of biomass (e.g., pine waste, forest wood residues, olive pomace) has its 
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Fig. 2 General schematic of a biomass SC 


own characterizing properties such as moisture content and heating value that 
determines its energy conversion efficiency. The pre-treatment block considers 
those activities that modify the quality (primarily moisture content) and/or shape 
of the biomass. Examples of this kind of processes are the chipping, pelletizing, 
drying, and torrefaction. Notice that these activities may be necessary, provided 
that there are some technologies which require feeding material to have a fixed 
maximum moisture content and/or some shape requirements in order to be 
processed. The generation block converts biomass into energy or any biofuel. 
Finally, the distribution block comprises those activities aiming at delivering the 
final product to the consumption points. 


2.1 Supply Chain Drivers 

All the aforementioned decisions will be taken such that an economic indicator, 
i.e., net present value (NPV), and an environmental impact metric, are optimised at 
the end of a predefined planning horizon. 


2.1.1 Cost Benefit Analysis 

The usage of NPV is the back bone of cost benefit analysis (CBA) of any project. 
A generic CBA consists of three steps: (i) valuation of the yearly costs and benefits 
of the project, (ii) discounting costs and benefits in future years to make them 
commensurate with present costs and benefits and (iii) calculation of the metric, 
(see Sect. 2.1 in chapter Global Clean Gas Process Synthesis and Optimisation). 
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Here, in order to compute the NPV, operational costs include those associated with 
production, handling of material, transportation and raw materials. Transportation 
costs are assumed to be linear functions of the actual flow of the product from the 
source echelon to the destination echelon. Revenue is obtained by the selling of 
products. Investments on facilities and equipment are also taken into account. 
The cash flows are discounted at a given return rate. 


2.1.2 Environmental Metrics 

Regarding the selection of environmental metrics, different methodologies have 
been developed; however, all of them rely on the accurate estimation of envi¬ 
ronmental interventions. Environment is compromised by industry mainly in two 
ways, namely, its emissions and the consumption of raw materials. See Sect. 2.3 in 
chapter Global Clean Gas Process Synthesis and Optimisation for further details. 

Here, the environmental metrics used are the ones devised in the work of 
Humbert et al. [17], which presents an implementation working at both mid-point 
and end-point (damage) levels. For each environmental intervention two charac¬ 
terisation factors are proposed, which eases model implementation. Their 
methodology, IMPACT 2002+, is mainly a combination between IMPACT 2002 
[18], Eco-indicator 99 [19] using egalitarian factors, CML [15], and the 
Intergovernmental Panel on Climate Change (IPCC) considerations for C0 2 
emissions. IMPACT 2002+ has grouped similar category end-points into a struc¬ 
tured set of damage categories by combining two main schools of impact model 
methods: classical impact assessment methods (CML/IPCC) and damage-oriented 
methods (Eco-indicator 99). This methodology proposes a feasible implementation 
of a combined mid-point/damage-oriented approach. It links all types of LCI results 
via 15 mid-point impacts (human toxicity, respiratory effects, ionizing radiation, 
ozone layer depletion, photochemical oxidation, aquatic ecotoxicity, terrestrial 
ecotoxicity, terrestrial acidification/nitrification, aquatic acidification, aquatic 
eutrophication, land occupation, global warming, non-renewable energy and min¬ 
eral extraction) to four areas of protection end-point categories (human health, 
ecosystem quality, climate change and global warming potential and resources). 

This approach contains the advantages of being able to calculate both mid- and 
end-point indicators. Within the presented model, and in order to avoid emission 
double counting, raw material emissions are not aggregated to product manufac¬ 
turing, similarly transport and energy consumption are considered separately. 

To conclude, Matthews et al. [20] highlight the importance of a footprint 
estimation that includes the total SC up to the production gate, also known as 
cradle-to-grave approach (tier 3). Furthermore, the authors refer to tier 4 emission 
estimations when the whole product life cycle is taken into account by considering 
emissions occurring during distribution and product end of life. This extended 
scope is expected to better aid effective environmental strategies since both firms 
and consumers have an important influence over the footprints through their 
‘purchase’ decisions. 
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3 The Mathematical Formulation 

This chapter describes, besides the methodology, also a tool that can be used to 
assist in the planning and design of a biomass SC under economical and 
environmental impact considerations. The resulting model is solved by using a 
multi-objective MILP (moMILP) algorithm, which allows observing possible 
environmental trade-offs between damage categories and the economic indicator. 

The mathematical formulation of the LCA-SC problem is briefly described 
next. The variables and constraints of the model can be roughly classified into 
three groups. The first one concerns process operations constraints given by the SC 
topology. The second one deals with the environmental model used while the third 
refers to the economic metric applied. 


3.1 Supply Chain Design: Planning Model 


The design-planning approach presented is a translation of the state task network 
(STN) formulation [21] to SC modelling, which is adapted from the work of 
Lamez et al. [22]. Such a formulation is suitable to collect all SC node information 
through a single variable, which eases the environmental and economic metrics 
formulation. This way SC node characteristics are modelled with a single equation 
set, since manufacturing nodes and distribution centres are treated in the same way 
as production and distribution activities. Subsequently, it turns out that the most 
important model variable is Pijg t \ which represents the specific activity of task 
i performed using technology j receiving input materials from site / and ‘deliv¬ 
ering’ output materials to site / during period t. Indeed, to model a production 
activity it must receive and deliver material within the same site (P^). In case of a 
distribution activity, facilities/and/ must be different. The model’s equations are 
briefly described in the next paragraphs. The separation between tasks and 
technologies allows for a flexible formulation of different scenarios. 

Materials mass balance must be satisfied at each one of the nodes. The expression 
for the mass balance for each material (state in the STN formulation) s consumed at 
each potential facility / in every time period t is presented next. Parameter oc si j is 
defined as the mass fraction of material s that is produced by task i performed using 
technology j; T s set refers to those tasks that have material s as output, while OL si j and 
T s set, refer to tasks that consume material s (Eq. 1). 

$sft—Sgft -1 = yy yy % syp^ft— yy yy v.?,/, t (i) 

/' i^Tsjefjipjf/) f y G (/.n/) 

The model assumes that process parameters are fixed (such as reaction 
conversion, separation factors and temperatures). This assumption is acceptable 
for most activities. Notice that the recipe for a given activity is fixed and given 
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by the parameters a si j and however, there are activities for which it may be 
desirable to let the model specify the mixture of inputs in order to achieve a 
given value for a specific biomass property (i.e., moisture content). For such 
activities the proportion of the different possible feedstocks should be variable. 
In order to account for those activities the mass balance shall be modified as 
shown in Eq. 2. 

Ssft ~ S s ft- 1 = CLsijPijfft — OLsijPijff't 

f ieT ‘jefrrjf,) f i£T s 

+ P y sijft — P v sijf't V^,/, t (2) 

ieT ° je(j t nj ff ) ief a je(JiPJf) 

For these flexible activities, it is necessary to make sure that the energy balance 
is achieved. This is done by introducing Eq. 3. Here, H \ si is the heating value for 
material s in activity i. Notice that heating values for feedstock are fixed. An 
activity changes the output heating value if (i) it is a pre-treatment task that 
modifies the biomass properties; or if (ii) it is a task that just changes the shape/ 
appearance of biomass but it is receiving different kinds of biomass as input. 

Y HV » Pv "ift = Y HV " Pv "jf< e 7 // t (3) 

s ^i seSi 

Let us consider that a flexible activity must accomplish a total moisture content. 
In such a case, constraint (4) must be satisfied. Parameters Water 5 and Water^ ax 
express the moisture content for material s and the maximum moisture content 
permitted for task i performed in equipment j. 

Y Water s Pv sijft < Water™* Y Pv sijft Mi e 7,7,/, t (4) 

seS i seSi 

Equation 5 models the temporal changes in facility capacities, in this sense the 
model allows for the simultaneous consideration of design and retrofit of SCs. 
Equation 6 serves for total capacity ( Fjf t ) bookkeeping taking into account the 
amount increased during planning period f(FE^) . 

VffiFEjfi < F% < VjftFEff MfJ € J f , t (5) 

Fjft = Fjf t -1 + FE^y-j, V /J 6 Jf, t (6) 

Equation 7 ensures the total production rate in each plant to be greater than a 
minimum desired production rate and lower than the available capacity. Further¬ 
more, parameter defines a minimum utilisation rate of technology j in site /, 
while Oijff determines the resource utilisation factor. 

PjfFjft-i — Oijff'Pfflt Ss Fjft- 1 

f ie/y 


V/,7 e J f ,t 


(7) 
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dijff' is the capacity utilisation rate of technology j by task i whose origin is 
location/and destination location f. This parameter is one of the key factors to be 
determined when addressing aggregated planning problems, considering strategic 
and tactical decisions. This operational model may be applied in continuous as 
well as in semi-continuous processes. First, let us consider the continuous pro¬ 
cesses. For these cases, the capacity utilisation factor is a conversion factor, 
which allows taking into account the equipment j capacity in site / in terms of 
task i kg of produced material per time unit. In this way, the 9ijff factor is the 
maximum throughput per planning period. On the other hand, this parameter is 
closely related to tasks operation time in the case of semi-continuous (batch) 
processes. Note that in this kind of production processes, the time period scale 
utilised in aggregated planning is usually larger than the time a task (production/ 
distribution activity) requires to be performed. Therefore, the sequencing-timing 
problem of short-term scheduling is transformed into a rough capacity problem 
where aggregated figures are used. It is important to have in mind that capacity 
is expressed as equipment j available time during one planning period, then 6^ 
represents the time required to perform task i in equipment j per unit of pro¬ 
duced material. Thus, once operation times are determined this parameter can be 
easily estimated. Equation 8 forces the amount of raw material s purchased from 
site / at each time period t to be lower than an upper bound given by physical 
limitations (A^). Also, the model assumes that part of the demand can actually 
be left unsatisfied because of limited production or supplier capacity. Thus, Eq. 9 
forces the sales of product s carried out in market / during time period t to be 


less than or equal to demand. 



Y, Y, Y Pffl't — ^sft 

Vs G RM,/ G Sup, t 

(8) 

f ieT s j£Ji 



Y Y Y Pyf'fi —^ en v 

Vs G FP,/ G Mkt, t 

(9) 


f i^Ts j^Ji 


For further model details the reader should refer to Lamez et al. [22]. 


3.2 Supply Chain: Environmental Model 

The application of the LCA methodology to a SC requires four steps, namely 
(i) goal setting, (ii) life cycle inventory (LCI), (iii) life cycle impact assessment 
(LCIA) and (iv) results interpretation towards improvement. 

Regarding goal setting, it is important to define the boundaries of the system 
under study, and which is the functional unit (FU) that the SC will provide. 
Boundaries in the case of the chemical industry in general and in the case of 
Biomass SCs are usually drawn from cradle-to-grave, this is due to the fact that 
most of these SC products (chemicals and electricity) are used in different ways 
and the use phase of products made from these products is too difficult to model 
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appropriately. Consequently raw material extraction, its processing and shipment 
to a market are considered as part of the SC system. Regarding the FU, commonly 
a certain amount of product produced is considered. In this sense it is advisable to 
compare different SCs in terms of the fulfilled amount of sales or portion of 
demand satisfied [12]. 

The LCI step requires the estimation of SC environmental interventions 
(emissions or natural raw material consumptions) which requires assessment of 
raw material producers, transportation and product processing impacts. This step is 
the most time consuming within a LCA due to the large amount of information that 
is required to be gathered, however the usage of LCI databases eases this issue. 
More importantly the use of a mathematical model helps in calculating the 
appropriate LCI for the optimal SC configuration. 

The results of the LCI step of the LCA can be interpreted by means of different 
environmental metrics. Environmental interventions are translated into metrics 
related to environmental impact as end-points or mid-points metrics by the usage 
of characterisation factors (CFs). This translation is the Life Cycle Impact 
Assessment (LCIA) step. The metrics used differ in their position along the 
environmental damage chain (environmental mechanism). 

The equations of the environmental model are briefly described next. 
Equation 10 models IC a f t which represents the mid-point a environmental impact 
associated to site / which rises from activities in period t\ is the a environ¬ 
mental category impact CF for task i performed using technology /, receiving 
materials from node/and delivering them at node/'. 

IC-aft = \ * ( 10 ) 

jeJf ieI j f 

Similarly to the case of a S y and a S y, the value of \is fixed and constant, 
provided that all environmental impacts are directly proportional to the activity 
performed in that node (Pijfft) . This issue is a common practice in LCA, where all 
direct environmental impacts are considered linear with respect to the FU [23]. In 
the case of transportation the FU commonly considered is the amount of material 
(kg) transported a given distance (kg km). Consequently the value of ican be 
calculated by Eq. 11 in the case of transportation, which considers the distance 
between sites (distance#*) and where i//Jj a represents the a environmental category 
impact CF for the transportation of a mass unit of material over a length unit. The 
study of environmental impacts associated to transport or production can be per¬ 
formed by setting the indices summation over the corresponding tasks (i.e., i E Tr 
or i E NTr). It should be noted that environmental impacts associated to materials 
transport are assigned to their origin node. 

> hjff'a = '/4 distance r V * G Tr j e J i - «>/>/' (! 1) 

Equation 12 introduces DamCg# which are a weighted sum of all mid-point 
environmental interventions combined using g end-point damage factors C ag and 
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then further normalised with NormFg factors. Equation 13 is used to compute the 
g normalised end-point damage along the whole SC ^DamC^. 

DamC, /( = Z NormF,,4 as ,IC a/( Vg,f,t (12) 

a£A g 

DamC® c = Z Z DamC^ Vg (13) 

f < 


Equations 14 and 15 sum up the environmental damage category results for 
each site / and for the whole SC, respectively. 

Impact 2002 = Z Z DamCg/j. V/ (14) 

8 * 

Impact 2 ™ 2 ,,, = Z Z Z DamC tf V/ (15) 

/ 8 * 


DamC^or Impact^^ii are both used as objective functions in the moMILP 
formulation. In this sense the use of damage categories is sometimes preferred 
given that they are easier to comprehend compared to mid-point values. 


3.3 Economic Model 

Many economic indicators have been proposed to assess the performance of a SC 
network design. The most traditional indicators are profit, NPV and total cost. Other 
more holistic measures have been recently proposed. Lamez et al. [24] proposed a 
model that pursues the maximisation of a financial key performance indicator, the 
corporate value of the firm at the end of the time horizon. The corporate value is 
computed by a discounted-free-cash-flow (DFCF) method which can be introduced 
as part of the mathematical formulation. Next, expressions to calculate (i) operating 
revenue, (ii) operating cost and (iii) capital investment are presented which would 
eventually permit integration with detailed financial models. Here, NPV will be 
used for the sake of simplicity and comprehensiveness. The application of other 
kind of metrics is out of the scope of this chapter, given the specific characteristics 
of the problem addressed in this work, but they are discussed under chapter 
Global Clean Gas Process Synthesis and Optimisation in Sect. 2.1 

Operating revenue is calculated by means of net sales which are the income 
source related to the normal SC activities. Thus, the total revenue incurred in any 
period t can be easily computed from products sales executed in period t as stated 
in Eq. 16. 


ESales^ = ^ ^ ^ Sales^Price^ \/t (16) 

sGFP/eMkt/'0(MktUSup) 
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In order to calculate overall operating cost an estimation of fixed costs and 
variable costs is required. The total fixed cost of operating a given SC structure can 
be computed using Eq. 17. Where FCFJ^ is the fixed unitary capacity cost of 
using technology j at site /. 

FCost ? = £ 5>CFJ ift F jft Vr (17) 

/^(MktUSup) jeJf 

The cost of purchases from supplier e , which is computed through Eq. 18, 
includes raw materials purchases, transport and production resources. 

EPurch e , = Purch™ + Purch* + Purch( ; rod Ve, t (18) 

The purchases (Purch™) associated to raw materials made to supplier e can be 
computed through Eq. 19. Variable / est represents the cost associated to raw 
material s purchased from supplier e. 

Purch™ = EE EE p mlest Ve € E rm , t (19) 

sEKMf £F e i£T s j^Ji 

The costs of transportation and production are determined by Eqs. 20 and 21, 
respectively. Here, p^, t denotes the e provider unitary transportation cost asso¬ 
ciated to material movement from location / to location f during period t. xJ^ et 
represents the unitary production cost associated to perform task i using tech¬ 
nology j, whereas x^ t represents the unitary inventory costs of material s storage at 
site /, both of them using provider e during period t. 


Purch *=E E 51V p W‘Pe/r't Vc € E tr , t 

ieTrjejiCJe f f 


( 20 ) 


purchr d = EE E E s ^t 

f i0r jeiJiHJf) s f${ SupUMkt) 


G Eprocb t 


( 21 ) 


In the case of x^ et , this parameter entails restrictions associated with oc si j and 
0 L si j, which forces the plant to operate at the same fixed conditions, meaning 
that the amount of utilities and labour spent is proportional to the amount of 
raw material processed. However the utilities and labour unitary cost may 
change over time. Moreover, possible cost decrease associated to economies of 
scale are disregarded by using the former assumption, higher production rates 
are associated linearly to higher production costs. Finally, the total investment 
on fixed assets is computed through Eq. 22. This equation includes the 
investment made to expand the technology’s capacity j in facility site / in 

period t (Price^FE 7 y^. 
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FAsset, = ^2^2 Price # F % + ( 22 ) 

/ J 

Equation 23 represents the calculation of profit at period t. To conclude, NPV is 
computed by means of Eq. 24. 


Profit^ = ESaleS; — 


FCost; + EPurch ej 


Vf 


(23) 



Profit, — FAsset, \ 
(1 + rate)' / 


(24) 


The selection of the discount rate (rate) for any time discounted metric is 
subject to controversy, given that it represents the trade-off between the enjoyment 
of present and future benefits and affects directly intergenerational aspects of 
sustainability. Higher values of rate devaluate future impacts and consequently 
they count little on long time horizon projects, which could be perceived as 
contrary to the interest of future generations. Identically to the case of a weighting 
set for a composite environmental index, the selection of a given discount rate is 
highly subjective and should represent the decision makers belief in terms of 
intergenerational aspects. 

Finally, the SC network design-planning problem whose objective is to opti¬ 
mise a given set of objective functions can be mathematically posed as follows: 

Min{-NPV, DamCf, Impact^} 


subject to Eqs. 2-24 

Xe {0,1};KM + 

Here X denotes the binary variables set, while Y corresponds to the continuous 
variable set. 


4 Case Study: A Biomass Supply Chain Geographically 
Located in Spain 

The case study used to illustrate the concepts behind the presented design strategy 
addresses a Biomass SC problem comparing the generation of electricity and H 2 
from two different kinds of feedstock: (i) different biomass wastes and (ii) coal, 
which is a wide extended fossil fuel resource. The SC under study comprises 
biomass collection sites and processing sites where biomass is pre-treated and 
used, in distribution centres and marketplaces. 
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Fig. 3 Location map for the potential SC network 


A simplified potential network is proposed and restricted to Spain (see Fig. 3). 
Lugo (FI), Ciudad Real (F2) and Burgos (F 3) are considered to be possible 
facilities location nodes. The feedstock is supposed to be available at Cordoba 
(LA), Lugo (LB), Cuenca (LC), Santander (LD) and Oviedo (LE). This last site is 
the one supplying coal. Hydrogen is supposed to be sold at three market places 
located at Madrid (Ml), Valencia (M2) and Barcelona (M3), while electricity is fed 
to the Spanish electricity network at their respective generation places. 

Different biomasses are modelled considering that each of them possesses 
different energy content and moisture. Table 1 shows the characteristics of the 
considered feedstocks, biomass wastes with useable energy content (see chapter 
Raw Materials, Selection, Preparation and Characterization) that can be found in 
Spain. Biomass types, approximate availability, costs and moisture contents are 
based on data published in Gomez et al. [25], Van Dyken et al. [10], Rentizelas 
et al. [8] and Panichelli and Gnansounou [4]. Lower heating values (LHV) for 
biomasses are taken from Phyllis database [26]. Hypothetically, coal has been 
considered as a dry material which does not require any pre-treatment. It is 
interesting to observe that higher moisture contents are present in forest and wood 
residues. Olive residues have low moisture content and high heating value which is 
comparable to the one of the sub-bituminous coal considered here. 

The biomass may be pre-treated before being finally processed. Figure 4 
depicts the different pre-treatment processes that may be applied to the biomass 
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Table 1 Feedstock properties. LHV is in ar (as received) basis; Wt: Weight 


Biomass 

Cost 
(€ ton - 

Lower 
x ) Heating 
Value, 
(MJ/kg) 

Moisture 

content 

(%wt) 

Seasonality 

Monthly 

availability 

(tons) 

Source 

location 

Forest wood residue 

25 

8.597 

48.9 

None 

40,500 

LA 

Pine waste 

35 

10.450 

40 

None 

3,350 

LB 

Almond tree prunings 

40 

11.313 

40 

December 
to February 

12,400 

LC 

Olive pomace 

65 

19.098 

7.6 

January to 
March 

73,400 

LD 

Olive pit 

35 

18.778 

6.1 

February to 
May 

27,900 

LD 

Coal 

45 

15.000 

- 

None 

65,000 

LE 


©-* 


Chipping 


.A 

-W BM-C I 



Pelletizer 

\ 



V 


/ 



Fig. 4 STN representing the pre-treatment activities for a generic biomass 

(BM) so that it achieves the adequate shape and properties (energy content and 
humidity) to follow up later processes. Several conditions have been assumed to 
select the specific path of pre-treatment units for each feedstock. They depend 
basically on the moisture content, shape and LHV of the input required for the 
treatment plant. Energy densification and matter densification through drying as 
well as dry matter loss in the different processes affecting the heating value have 
been considered. Nevertheless, in this case study BM bulk density has not been 
considered for the sake of simplicity and also given the planning horizon required 
for the strategic decisions that are being addressed. Note that coal does not need 
any pre-treatment. The characteristics for the pre-treatment processes are pre¬ 
sented in Table 2 The data are from Panichelli and Gnansounou [4] and from 
Hamelinck et al. [5]. The pre-treatment options considered here are: 

• Chipper. Transformation of the biomass as received, into chips. It is the first 
biomass size reduction step, and mandatory for all the biomass sources con¬ 
sidered. It is considered that a first step of shape homogenisation is crucial for an 
integrated Biomass SC. Moisture content is not modified, but there exists a loss 
of dry matter, which is considered. 
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Table 2 Pre-treatment processes and their main modelling assumptions 


Activity/ 

equipment 

Moisture 
losses (%) 

Dry matter 
losses (%) 

Operating 
cost (€ ton - 

Capacity 
l ) (tons h -1 ) 

Investment 
(1 x 10 6 €) 

Electricity 
consumption 
(MWh ton -1 ) 

Chipper 

0 

0.17 

2.5 

30 

0.370 

5 

Dryer 

88 

0.08 

55 

100 

5.000 

20 

Torref actor 

55 

19 

40 

20 

0.100 

37 

Pelletizer 

0 

0 

3.5 

6 

0.485 

30 


Table 3 Parameters for the processes for electricity and H 2 

generation 


Technology 

Operating 
cost (€) 

Capacity 

Investment 
(1 x 10 6 €) 

Product price 
(€) 

Total monthly 
demand 

Electricity 

34.2 MWh -1 

300 MW 

860 

0.151 kWh -1 

75,000 MWh 

h 2 

1,880 ton' 1 

33.6 tons h -] 

' 1,500 

3 kg- 1 

650 ton 


• Dryer. Active dryer is needed when the source has too high moisture contents. The 
condition to enter this step is to have a humidity higher than 7%. Passive drying is 
assumed not to be significant according to the unit of time that the model con¬ 
siders. The dry matter loss is insignificant if compared with the moisture loss. 

• Torrefactor. Its main objective is the LHV increase through volatiles release in 
an inert atmosphere. The biggest dry matter loss is obtained here, and the drying 
process also takes place. The condition to pass through this unit is to have a 
LHV lesser than 15 MJ/kg. 

• Pelletizer. In this case study we assume that the biomass must be pelletized 
before its final transformation, thus a high homogeneous raw material is arriving 
to the final plants. Note that the pelletizer requires input to have humidity equal 
to or lower than 7%. This requirement may be fulfilled by mixing different types 
of biomasses or pre-treated material. It is the main characteristic of this unit 
(note that no moisture or dry matter losses are considered here). 

The technology that is employed to provide the final product is gasification. An 
integrated gasification combined cycle (IGCC) power plant is assumed for the 
electricity generation, and a gasification plant with carbon capture and storage 
(CCS) for the H 2 generation. Efficiencies of 40 and 30% are assumed for each plant, 
respectively, the adoption of such figures will be extensively discussed in chapter 
Modelling Superstructure for Conceptual Design of Syngas Generation and 
Treatment. It is assumed that the lower capacity that can be installed for the 
energy generator is 85 MWh. Other relevant information concerning these tech¬ 
nologies is presented in Table 3. These data are from Perez-Fortes et al. [27] and 
IEA-GHG [28]. Biomass and H 2 transportation prices were estimated at 0.03 and 
0.05 € ton -1 km -1 from current economical trends, respectively. The demand of H 2 
is evenly distributed along the three markets (Ml, M2 and M3, while the demand of 
electricity is supplied to the grid from any facility location). It is assumed that the 
demand must be completely satisfied. For the case of NPV optimisation return rate is 
assumed to be 8% which is a common value for this kind of projects. 








44 


J. M. Larnez et al. 


In order to assess the environmental impact associated to the biomass-based 
energy SC, the available LCI values were retrieved from the LCI database 
EcoinventV1.3 [29] and using SimaPro 7.1.6 [30] and converted directly to the 
IMPACT 2002+ mid-point indicators. For those activities which were not available, 
the impacts were assumed based on similar products or activities. The environ¬ 
mental impacts associated to energy generation, H 2 production and pre-treatment 
processes without consideration of feedstock and transportation consumption are 
found in Table 4. The environmental impact for transportation activities is pre¬ 
sented in that table as well. The environmental impact for feedstock can be found in 
Table 5 which does not consider impacts associated to transportation. 

The project is evaluated along a planning horizon of 25 years, considering 
monthly planning decisions. The model has been implemented in GAMS which is 
an algebraic modelling software. The formulation of the SC-LCA model leads to a 
MILP with 4,159 equations, 41,221 continuous variables and 96 discrete variables. 
It takes 61 CPU seconds to reach a solution with a 0.1% integrality gap on a 
2.0 GHz Intel Core 2 Duo computer using the MIP solver of CPLEX. 

Figure 5 shows the obtained dominant biomass-based SC that maximises NPV. 
It is found that the three potential locations are considered and on each one of them 
a facility is opened. All pre-treatment technologies are installed in location FI 
besides the required equipment to produce H 2 . From this site H 2 is delivered to all 
markets. Note that FI is collecting all the forest wood residues (FWR) for which 
larger mass flows are required due to their low FHV. By establishing FI, which is 
near to the FWR collection site, significant savings in transportation are obtained. 
The electricity is generated in site F2. In this site; equipments to perform chipping, 
drying and pelletising are installed. The electricity demand of each market is 
satisfied from site F2. Site F3 is used just as a distribution centre for pre-treated 
biomass. Equipment for chipping and drying is installed in such a site. 

Table 6 shows the proposed capacity to be installed for each of the equipments 
at every site to obtain the maximum NPV configuration. Notice that for this 
configuration there are some inter-site flows, clearly showing the capabilities of 
the model to tackle with inter-site distribution tasks. Forest wood residues which 
have been dried and torrefied are being sent from site FI to F2 , while F3 is 
transferring dried pine waste and dried almond tree pruning to location F2 in order 
to be converted to energy later on. By having material flows of pre-treated biomass 
the transportation cost is reduced due to the mass decrease that is achieved through 
the utilisation of such processes. 

The optimal configuration for the environmental impact has also been obtained. 
Figure 6 shows the minimum IMPACT 2002+ configuration for the biomass-based 
SC. Please note that this supply chain fulfills with the same demand as the one 
obtained by optimising NPV. The capacity proposed to be installed for the 
equipment of this configuration is presented in Table 7. Note that for this case 
the location F3 is not considered, and all biomass is sent from the collection sites 
to locations FI and F2. This configuration is satisfying the demand of electricity 
from both locations FI and F2 , whereas H 2 is delivered from site F2. There is one 
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Fig. 5 Optimum NPV 
network configuration for the 
biomass-based SC 



Table 6 Equipment capacity 
for the optimum NPV 
Biomass SC 


Equipment 

Facility 



FI 

F2 

F3 

Chipper (tons h -1 ) 

67.5 

60.4 

40.0 

Dryer (tons h -1 ) 

47.4 

20.7 

20.0 

Torrefactor (tons h _1 ) 

55.8 

- 

- 

Pelletizer (tons h -1 ) 

23.2 

60.3 

- 

Electricity (GJ h _1 ) 


450.0 

- 

H 2 (tons h -1 ) 

1.1 


- 


inter-site flow from site F2 to FI , which corresponds to a flow of dried olive 
pomace, in this environmental friendly configuration. 

Recall that we introduced a ‘flexible’ task to account for those tasks for which 
we would like the model to decide how to better mix different biomasses so as to 
achieve a given specified biomass property. We have assumed that the pelletizer is 
one of such tasks for this case study. To give an example, there are periods in 
which the model proposes to make the following mix: 1.4% forest wood residues, 
30.3% dried and torrefied forest wood residues, 10.5% dried pine waste, 14.4% 
dried almond tree prunings and 43.5% chipped olive pomace (mass basis). This 
mixture is then fed to the syngas production plant. The values of humidity cor¬ 
responding to these materials are 10.0, 6.0, 7.0, 7.0 and 7.5%, respectively. It can 
be proved that the humidity of this mix is 7.0% which is the maximum humidity 
allowed for the pelletizer input. 
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Fig. 6 Optimum Impact 
2002+ network configuration 
for the biomass-based SC 



Table 7 Equipment capacity 
for the optimum Impact 
2002+ Biomass SC 


Equipment 

Facility 



FI 

F2 

F3 

Chipper (tons h -1 ) 

73.1 

70.3 

- 

Dryer (tons h _ ) 

54.4 

20.7 

- 

Torrefactor (tons h _1 ) 

51.7 

- 

- 

Pelletizer (tons h -1 ) 

40.0 

62.7 

- 

Electricity (GJ h _1 ) 

300.0 

300.0 


H 2 (tons h -1 ) 


1.1 



Tables 8 and 9 summarize the most significant values corresponding to both 
solutions regarding environmental and economic aspects. By deploying the SC 
configuration corresponding to the more profitable SC configuration, a NPV equal 
to 228.51 M€ is obtained. This value is reduced by 3% when the environmental 
friendly configuration is established. Note that the main difference between these 
two configurations is the investment required for installing the proposed capacity 
in the different sites. For the case of NPV maximisation the installation of an 
electricity generation plant at location F2 makes it necessary to send larger amount 
of biomass to this site. In order to reduce this additional cost a distribution centre 
at location F3 is established where biomass is treated before being transferred to 
the other sites. The reduction in capacity investment offsets, the increase in the 
transportation cost and the extra costs are associated to the pre-treatment of bio¬ 
mass. Table 8 also shows the payback periods and internal rate of return for both 
solutions. 
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Table 8 Economic aspects 
arising from single objective 

Impact 2002+ NPV 

optimisation optimisation 

upuim^ciLiuii v diiu 




Impact 2002+) (M€) 

Investment 

170.21 

148.13 


Biomass cost 

820.16 

819.87 


Transportation cost 

92.73 

122.65 


Production cost 

2,222.77 

2,224.02 


Sales 

3,987.00 

3,987.00 


Profit 

851.34 

820.46 


NPV a 

220.59 

228.51 


IRR (%) b 

24.02 

27.38 


Payback period (years) 

5.00 

4.51 


a Based on a return rate of 8% 



b Based on a planning horizon of 25 years 


Table 9 Environmental 

End-point impact 

Impact 2002+ 

NPV optimisation 

impacts arising from single 

category 

optimisation 


economic and overall 
environmental objective 

Human health 

16,255.29 

17,267.21 

function optimisation results Ecosystem quality 

3,375.79 

3,610.96 

(Impact 2002+ points per 

Climate change 

90,383.37 

90,950.66 

year) 

Resources 

5,292.64 

5,852.73 


Impact 2002+ 

115,307.09 

117,681.56 


Fig. 7 Distribution of annual 
environmental impacts for 
single objective optimisation 
solutions, according to 
different SC activities 



Impact 2002+ minimization NPV maximization 


■ Production 

■ Pretreatment 

■ Product transportation 

■ Raw Biom. transportation 

■ Biomass 


With regard to environmental interventions, the NPV optimum solution renders 
an environmental impact of 117681 pts. The environmental friendly configuration 
can slightly reduce the impact by 2%. Notice that the impacts for each damage 
category (see Table 9) are very similar for both solutions. These configurations 
have a major impact on the climate change category which represents 78% of the 
overall impact. Figure 7 presents the distribution of environmental impacts 
according to different SC activities. Notice that as for the NPV comparison, the 
difference between the environmental impacts is mainly due to transportation. 
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Fig. 8 Energy price versus 0155 “ 
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In the presented case study electricity generation and H 2 production are the most 
important factors contributing to the overall environmental impact in both single 
objective optimisation cases; while biomass sourcing is the least impacting aspect. 
This clearly shows that activities to reduce environmental impact should be 
focused on improving the technologies used to produce energy and H 2 . 

We would like to highlight how sensitive the solutions for the biomass-based 
SC are to the final product prices. For instance if we assume that the energy price 
is reduced to 50% (i.e., 0.075 € kWh -1 ) the NPV at the end of the planning 
horizon (i.e., 25 years) would be negative and equal to —553.25 x 10 6 €. The 
price that renders an internal rate of return of 8% considering a project life time of 
25 years is 0.129 € kWh -1 . Any price below this value would require some sort of 
subsidy. Figure 8 shows the energy prices that are required for different internal 
rates of return. 

For comparison purposes the optimal SC based on coal has been also obtained. 
The optimal NPV configuration for this case proposes to deliver electricity from 
location FI , where a capacity of 450 GJ h -1 is installed, while the H 2 is produced 
at location F3 where a capacity of 1.1 tons h -1 is installed. The minimum Impact 
2002+ solution for this case (where coal is the only raw material available) is the 
same as the one obtained by the NPV optimisation. The reason is that the only way 
to improve the environmental impacts is by means of transportation, which is the 
case of NPV optimisation since there are no pre-treatment associated with coal. 
This solution is summarised in Tables 10 and 11. Note that an NPV improvement 
of 219% can be gained by utilising coal as feedstock when compared to the 
biomass-based SC. Notice that the main difference is from the production cost 
which is due to the pre-treatment activities that are required in the biomass-based 
SC. This fact also makes the investment increase in the Biomass SC. 

Regarding the environmental impacts, the Impact 2002+ is increased in 203% 
compared with the biomass-based SC. It is noteworthy that the impact associated 
with the climate change category is very similar for both cases. We have to bear in 
mind that C0 2 is still emitted when using a Biomass SC, however this biomass is 
regenerated faster than fossil fuels. Nevertheless, the other categories are signifi¬ 
cantly increased specially for the case of resources and human health. The resource 
regeneration issue is the reason why the impact related to the resources category is 
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Table 10 Economic aspects 
arising from single objective 


NPV optimisation 

Investment 

111.74 

optimisation for the coal- 

based SC (M€) 

Coal cost 

839.55 


Transportation cost 

181.70 


Production cost 

1,136.10 


Sales 

3,987.00 


Profit 

1,829.65 


NPV a 

729.40 


IRR (%) b 

112.75 


Payback period (years) 

1.53 

Table 11 Environmental 
impacts arising from single 

a Based on a return rate of 8% 
b Based on a planning horizon of 25 years 


End-point impact category 

NPV optimisation 

Human health 

109,640.53 

economic and overall 

environmental objective 

Ecosystem quality 

11,077.16 

function optimisation results 

Climate change 

95,334.48 

for the coal-based SC (Impact 

Resources 

140,605.89 

2002+ points per year) 

Impact 2002+ 

356,658.06 


Fig. 9 Distribution of annual 
environmental impacts for the 
coal-based SC 
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NPV optimization 


■ Product 
transportation 

□ Coal transportation 

■ Coal sourcing 

□ Production 


larger in the coal-based case. This fact also emphasises the importance of having 
an overall impact indicator instead of a partial indicator such as C0 2 kg. Notice 
that if we compare the coal and biomass-based SCs based on C0 2 kg, there would 
not be an important reduction in the emissions to environmentally justify the 
project. Figure 9 shows the distribution of environmental impacts according to 
different SC activities for the coal-based SC. Notice that for the coal-based case 
the impact of production accounts for 46% of the overall impact. The other activity 
generating environmental impacts is the coal sourcing which accounts for the 50% 
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of the overall impact. Moreover, from Fig. 9, it can also be concluded that if the 
sourcing of raw materials would be disregarded then different solutions would be 
obtained. If we disregard the raw materials sourcing from the analysis, the bio- 
mass-based SC and the coal-based SC would result in a environmental impact of 
112,701 and 178,225 points, respectively. This highlights the significance of a SC 
approach to the problem, rather than focusing only in the processing sites. This 
also points out the relevance of ‘purchase’ decisions on the environmental impact 
of an SC. 


5 Conclusions 

This chapter presents an approach for designing and planning efficient Biomass 
SCs. The model consists of a multi-period MILP that accounts for the multi¬ 
objective optimisation of economic and environmental interventions. The model 
considered the long-term strategic decisions (e.g., establishing of pre-treatment 
trains of units and their respective location, selection of biomass sources, location 
of processing sites and distribution centres). The problem has been formulated as a 
multi-objective optimisation, where two objective functions are considered, the net 
present value and the Impact 2002+ metric was adopted as a measure of overall 
environmental impact. 

A biomass SC case study, geographically located in Spain, is presented here 
considering the availability of five different potential biomass waste sources. It has 
been considered that two final products, electricity and H 2 , are delivered to mar¬ 
ketplaces by means of gasification plants. The most environmental friendly and the 
most profitable SC configurations for this case study have been shown and their 
differences have been discussed. However, the indicator values were very similar. It 
was shown that the more impacting damage category was climate change which 
accounts for approximately 80% of the overall environmental impact. Moreover, the 
sensitiveness of the optimal solutions to change in prices and rate of return values 
was demonstrated. The most profitable SC results in a positive NPV, however 
reductions of product prices can easily lead to economic losses. The prices to achieve 
different internal rates over a project life time of 25 years were presented. For this 
specific case study it was demonstrated that if a return rate of 8%, which is typical for 
this type of projects, is required, the price cannot be lower than 0.129 € kWh -1 . 
Otherwise the project should be subject to a subsidy if it is meant to be sustainable. 
This demonstrates how this type of models can be used to determine subsidy policies 
in order to actually drive industry towards more environmental practices. 

The problem was also solved for a coal-based SC. This comparison empha¬ 
sised the relevance of (i) using an overall environmental impact rather than a 
partial metric such as kg of C0 2 , and (ii) a supply chain approach to the problem 
so that the entire life cycle of the product and/or process is included and ana¬ 
lysed. Narrower approaches (i.e., merely analyzing the processing/manufacturing 
sites) may lead to biased solutions when the magnitude of environmental impacts 
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associated with producing, collecting or extracting the raw materials are con¬ 
siderably high. 
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